﻿<?xml version="1.0" encoding="utf-8" standalone="yes"?>
<?application progid="SMath Studio Desktop" version="0.98.6179.21440"?>
<regions>
  <settings>
    <identity>
      <id>09482963-2a40-416e-a67b-1c04f466f7c6</id>
      <revision>19</revision>
    </identity>
    <calculation>
      <precision>3</precision>
      <exponentialThreshold>5</exponentialThreshold>
      <fractions>auto</fractions>
    </calculation>
    <pageModel active="true" printAreas="true" simpleEqualsOnly="false" printBackgroundImages="true">
      <paper id="1" orientation="Portrait" width="850" height="1100" />
      <margins left="39" right="39" top="39" bottom="39" />
      <header alignment="Center" color="#a9a9a9">&amp;[DATE] &amp;[TIME] - &amp;[FILENAME]</header>
      <footer alignment="Center" color="#a9a9a9">&amp;[PAGENUM] / &amp;[COUNT]</footer>
      <backgrounds />
    </pageModel>
    <dependencies>
      <assembly name="SMath Studio Desktop" version="0.98.6179.21440" guid="a37cba83-b69c-4c71-9992-55ff666763bd" />
      <assembly name="Math Region" version="0.98.6179.21440" guid="02f1ab51-215b-466e-a74d-5d8b1cf85e8d" />
      <assembly name="Special Functions" version="1.11.6179.21442" guid="2814e667-4e12-48b1-8d51-194e480eabc5" />
      <assembly name="Text Region" version="1.10.6179.21446" guid="485d28c5-349a-48b6-93be-12a35a1c1e39" />
      <assembly name="Plot Region" version="1.9.6179.21450" guid="c451c2b5-798b-4f08-b9ec-b90963d1ddaa" />
    </dependencies>
  </settings>
  <region id="0" top="18" color="#ff0000" bgColor="#ffffff">
    <area collapsed="true">
      <title lang="eng">
        <p>     XY data set     </p>
      </title>
    </area>
    <region id="1" left="9" top="54" width="110" height="567" color="#000000" bgColor="#ffffff" fontSize="10">
      <math>
        <input>
          <e type="operand">XY</e>
          <e type="operand">1</e>
          <e type="operand">4.18</e>
          <e type="operand">2</e>
          <e type="operand">4.67</e>
          <e type="operand">3</e>
          <e type="operand">5.3</e>
          <e type="operand">4</e>
          <e type="operand">5.37</e>
          <e type="operand">5</e>
          <e type="operand">5.45</e>
          <e type="operand">6</e>
          <e type="operand">5.74</e>
          <e type="operand">7</e>
          <e type="operand">5.65</e>
          <e type="operand">8</e>
          <e type="operand">5.84</e>
          <e type="operand">9</e>
          <e type="operand">6.36</e>
          <e type="operand">10</e>
          <e type="operand">6.38</e>
          <e type="operand">11</e>
          <e type="operand">6.31</e>
          <e type="operand">12</e>
          <e type="operand">6.46</e>
          <e type="operand">13</e>
          <e type="operand">6.48</e>
          <e type="operand">14</e>
          <e type="operand">6.59</e>
          <e type="operand">15</e>
          <e type="operand">6.74</e>
          <e type="operand">16</e>
          <e type="operand">6.56</e>
          <e type="operand">17</e>
          <e type="operand">6.85</e>
          <e type="operand">18</e>
          <e type="operand">6.53</e>
          <e type="operand">19</e>
          <e type="operand">6.68</e>
          <e type="operand">20</e>
          <e type="operand">7.01</e>
          <e type="operand">21</e>
          <e type="operand">6.82</e>
          <e type="operand">22</e>
          <e type="operand">6.65</e>
          <e type="operand">23</e>
          <e type="operand">6.95</e>
          <e type="operand">24</e>
          <e type="operand">7.12</e>
          <e type="operand">25</e>
          <e type="operand">7.14</e>
          <e type="operand">26</e>
          <e type="operand">6.91</e>
          <e type="operand">27</e>
          <e type="operand">6.94</e>
          <e type="operand">28</e>
          <e type="operand">7.23</e>
          <e type="operand">29</e>
          <e type="operand">6.9</e>
          <e type="operand">30</e>
          <e type="operand">7.39</e>
          <e type="operand">31</e>
          <e type="operand">6.94</e>
          <e type="operand">31</e>
          <e type="operand">2</e>
          <e type="function" preserve="true" args="64">mat</e>
          <e type="operator" args="2">:</e>
        </input>
      </math>
    </region>
    <region id="2" top="612" color="#000000" bgColor="#ffffff">
      <area terminator="true" />
    </region>
  </region>
  <region id="3" top="648" color="#ff0000" bgColor="#ffffff">
    <area collapsed="true">
      <title lang="eng">
        <p>     linfit modules     </p>
      </title>
    </area>
    <region id="4" left="9" top="684" width="595" height="142" color="#000000" bgColor="#ebebeb" fontSize="10">
      <math fractionType="none" decimalPlaces="4" exponentialThreshold="3">
        <input>
          <e type="operand">vx</e>
          <e type="operand">vy</e>
          <e type="operand">a</e>
          <e type="operand">b</e>
          <e type="operand">c</e>
          <e type="function" args="5">plotG</e>
          <e type="operand">a</e>
          <e type="operand" style="string">char</e>
          <e type="operator" args="2">≡</e>
          <e type="operand">b</e>
          <e type="operand">size</e>
          <e type="operator" args="2">≡</e>
          <e type="operand">c</e>
          <e type="operand" style="string">color</e>
          <e type="operator" args="2">≡</e>
          <e type="operand">1</e>
          <e type="operand">3</e>
          <e type="function" preserve="true" args="5">mat</e>
          <e type="operand">n</e>
          <e type="operand">vx</e>
          <e type="function" preserve="true" args="1">length</e>
          <e type="operator" args="2">:</e>
          <e type="operand">plot</e>
          <e type="operand">vx</e>
          <e type="operand">1</e>
          <e type="function" preserve="true" args="2">el</e>
          <e type="operand">vy</e>
          <e type="operand">1</e>
          <e type="function" preserve="true" args="2">el</e>
          <e type="operand">a</e>
          <e type="operand">b</e>
          <e type="operand">c</e>
          <e type="function" preserve="true" args="5">augment</e>
          <e type="operator" args="2">:</e>
          <e type="operand">i</e>
          <e type="operand">2</e>
          <e type="operand">n</e>
          <e type="function" preserve="true" args="2">range</e>
          <e type="operand">plot</e>
          <e type="operand">plot</e>
          <e type="operand">vx</e>
          <e type="operand">i</e>
          <e type="function" preserve="true" args="2">el</e>
          <e type="operand">vy</e>
          <e type="operand">i</e>
          <e type="function" preserve="true" args="2">el</e>
          <e type="operand">a</e>
          <e type="operand">b</e>
          <e type="operand">c</e>
          <e type="function" preserve="true" args="5">augment</e>
          <e type="function" preserve="true" args="2">stack</e>
          <e type="operator" args="2">:</e>
          <e type="function" preserve="true" args="3">for</e>
          <e type="operand">plot</e>
          <e type="operand">5</e>
          <e type="operand">1</e>
          <e type="function" preserve="true" args="7">line</e>
          <e type="operator" args="2">:</e>
        </input>
      </math>
    </region>
    <region id="5" left="9" top="828" width="435" height="365" border="true" color="#000000" bgColor="#ebebeb" fontSize="10">
      <math optimize="2" fractionType="fraction" decimalPlaces="4">
        <description active="true" position="Top" lang="eng">
          <p>Applicable parts are: Optimiz "numeric"</p>
        </description>
        <input>
          <e type="operand">XY</e>
          <e type="operand">f</e>
          <e type="operand">φ</e>
          <e type="operand">β</e>
          <e type="function" args="4">βlinfit</e>
          <e type="operand">X</e>
          <e type="operand">XY</e>
          <e type="operand">1</e>
          <e type="function" preserve="true" args="2">col</e>
          <e type="operator" args="2">:</e>
          <e type="operand">Y</e>
          <e type="operand">XY</e>
          <e type="operand">2</e>
          <e type="function" preserve="true" args="2">col</e>
          <e type="operator" args="2">:</e>
          <e type="operand">n</e>
          <e type="operand">XY</e>
          <e type="function" preserve="true" args="1">rows</e>
          <e type="operator" args="2">:</e>
          <e type="operand">k</e>
          <e type="operand">x</e>
          <e type="operand">β</e>
          <e type="function" args="2">φ</e>
          <e type="function" preserve="true" args="1">length</e>
          <e type="operand">1</e>
          <e type="operator" args="2">-</e>
          <e type="operator" args="2">:</e>
          <e type="operand" style="string"> vx raw initial vector fit </e>
          <e type="operand">i</e>
          <e type="operand">1</e>
          <e type="operand">X</e>
          <e type="function" preserve="true" args="1">length</e>
          <e type="function" preserve="true" args="2">range</e>
          <e type="operand">vx</e>
          <e type="operand">i</e>
          <e type="function" preserve="true" args="2">el</e>
          <e type="operand">X</e>
          <e type="operand">i</e>
          <e type="function" preserve="true" args="2">el</e>
          <e type="operand">β</e>
          <e type="function" args="2">f</e>
          <e type="function" preserve="true" args="1">eval</e>
          <e type="operator" args="2">:</e>
          <e type="function" preserve="true" args="3">for</e>
          <e type="operand" style="string">M-&gt; Cholesky Gradient Mat</e>
          <e type="operand">j</e>
          <e type="operand">1</e>
          <e type="operand">k</e>
          <e type="function" preserve="true" args="2">range</e>
          <e type="operand">i</e>
          <e type="operand">1</e>
          <e type="operand">n</e>
          <e type="function" preserve="true" args="2">range</e>
          <e type="operand">M</e>
          <e type="operand">i</e>
          <e type="operand">j</e>
          <e type="function" preserve="true" args="3">el</e>
          <e type="operand">X</e>
          <e type="operand">i</e>
          <e type="function" preserve="true" args="2">el</e>
          <e type="operand">β</e>
          <e type="function" args="2">φ</e>
          <e type="operand">j</e>
          <e type="operand">1</e>
          <e type="operator" args="2">+</e>
          <e type="function" preserve="true" args="2">el</e>
          <e type="function" preserve="true" args="1">eval</e>
          <e type="operator" args="2">:</e>
          <e type="function" preserve="true" args="3">for</e>
          <e type="function" preserve="true" args="3">for</e>
          <e type="operand" style="string">Cholesky solver [LS]</e>
          <e type="operand">Φ</e>
          <e type="operand">M</e>
          <e type="function" preserve="true" args="1">transpose</e>
          <e type="operand">M</e>
          <e type="operator" args="2">*</e>
          <e type="operator" args="2">:</e>
          <e type="operand">d</e>
          <e type="operand">M</e>
          <e type="function" preserve="true" args="1">transpose</e>
          <e type="operand">Y</e>
          <e type="operand">vx</e>
          <e type="operator" args="2">-</e>
          <e type="bracket">(</e>
          <e type="operator" args="2">*</e>
          <e type="operator" args="2">:</e>
          <e type="operand">Δβ</e>
          <e type="operand">Φ</e>
          <e type="operand">1</e>
          <e type="operator" args="1">-</e>
          <e type="operator" args="2">^</e>
          <e type="function" preserve="true" args="1">eval</e>
          <e type="operand">d</e>
          <e type="operator" args="2">*</e>
          <e type="operator" args="2">:</e>
          <e type="operand">b</e>
          <e type="operand">β</e>
          <e type="operand">Δβ</e>
          <e type="operator" args="2">+</e>
          <e type="operator" args="2">:</e>
          <e type="operand">13</e>
          <e type="operand">1</e>
          <e type="function" preserve="true" args="15">line</e>
          <e type="operator" args="2">:</e>
        </input>
      </math>
    </region>
    <region id="6" left="441" top="828" width="330" height="201" border="true" color="#000000" bgColor="#ebebeb" fontSize="10">
      <math fractionType="fraction" decimalPlaces="4">
        <description active="true" position="Top" lang="eng">
          <p>Vector gradient ... Optimiz Symbolic </p>
        </description>
        <input>
          <e type="operand">x</e>
          <e type="operand">β</e>
          <e type="operand">f</e>
          <e type="operand">n</e>
          <e type="function" args="4">φ</e>
          <e type="operand">ct</e>
          <e type="operand">1</e>
          <e type="operator" args="2">:</e>
          <e type="operand">i</e>
          <e type="operand">1</e>
          <e type="operand">n</e>
          <e type="function" preserve="true" args="2">range</e>
          <e type="operand">PD</e>
          <e type="operand">ct</e>
          <e type="function" preserve="true" args="2">el</e>
          <e type="operand">x</e>
          <e type="operand">β</e>
          <e type="function" args="2">f</e>
          <e type="bracket">(</e>
          <e type="operand">β</e>
          <e type="operand">ct</e>
          <e type="function" preserve="true" args="2">el</e>
          <e type="function" preserve="true" args="2">diff</e>
          <e type="operator" args="2">:</e>
          <e type="operand" style="string">stack f(x,β) &amp; PD </e>
          <e type="operand">V</e>
          <e type="operand">1</e>
          <e type="function" preserve="true" args="2">el</e>
          <e type="operand">x</e>
          <e type="operand">β</e>
          <e type="function" args="2">f</e>
          <e type="operator" args="2">:</e>
          <e type="operand">V</e>
          <e type="operand">i</e>
          <e type="operand">1</e>
          <e type="operator" args="2">+</e>
          <e type="function" preserve="true" args="2">el</e>
          <e type="operand">PD</e>
          <e type="operand">i</e>
          <e type="function" preserve="true" args="2">el</e>
          <e type="operator" args="2">:</e>
          <e type="operand">2</e>
          <e type="operand">1</e>
          <e type="function" preserve="true" args="4">line</e>
          <e type="operand">ct</e>
          <e type="operand">ct</e>
          <e type="operand">1</e>
          <e type="operator" args="2">+</e>
          <e type="operator" args="2">:</e>
          <e type="operand">4</e>
          <e type="operand">1</e>
          <e type="function" preserve="true" args="6">line</e>
          <e type="function" preserve="true" args="3">for</e>
          <e type="operand">V</e>
          <e type="operand">3</e>
          <e type="operand">1</e>
          <e type="function" preserve="true" args="5">line</e>
          <e type="operator" args="2">:</e>
        </input>
      </math>
    </region>
    <region id="7" top="1233" color="#000000" bgColor="#ffffff">
      <area terminator="true" />
    </region>
  </region>
  <region id="8" left="36" top="1269" width="395" height="31" color="#000000" bgColor="#80ffff" fontSize="14">
    <text lang="eng">
      <p bold="true">Part 1: Dress the "linfit" solver</p>
    </text>
  </region>
  <region id="9" left="36" top="1305" width="295" height="69" color="#000000" bgColor="#ffffff" fontSize="10">
    <math>
      <input>
        <e type="operand">X</e>
        <e type="operand">XY</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="2">col</e>
        <e type="operator" args="2">:</e>
        <e type="operand">Y</e>
        <e type="operand">XY</e>
        <e type="operand">2</e>
        <e type="function" preserve="true" args="2">col</e>
        <e type="operator" args="2">:</e>
        <e type="operand">data</e>
        <e type="operand">X</e>
        <e type="operand">Y</e>
        <e type="operand" style="string">o</e>
        <e type="operand">5</e>
        <e type="operand" style="string">black</e>
        <e type="function" args="5">plotG</e>
        <e type="operator" args="2">:</e>
        <e type="operand">3</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="5">line</e>
      </input>
    </math>
  </region>
  <region id="10" left="36" top="1395" width="154" height="88" color="#000000" bgColor="#ffffff" fontSize="10">
    <math>
      <description active="true" position="Top" lang="eng">
        <p>1. model combination</p>
      </description>
      <input>
        <e type="operand">x</e>
        <e type="operand">β</e>
        <e type="function" args="2">p</e>
        <e type="operand">β</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="2">el</e>
        <e type="operand">β</e>
        <e type="operand">2</e>
        <e type="function" preserve="true" args="2">el</e>
        <e type="operand">x</e>
        <e type="function" preserve="true" args="1">ln</e>
        <e type="operator" args="2">*</e>
        <e type="operand">β</e>
        <e type="operand">3</e>
        <e type="function" preserve="true" args="2">el</e>
        <e type="operand">x</e>
        <e type="function" preserve="true" args="1">sqrt</e>
        <e type="operator" args="2">*</e>
        <e type="operand">3</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="5">mat</e>
        <e type="operator" args="2">:</e>
      </input>
    </math>
  </region>
  <region id="11" left="216" top="1422" width="192" height="24" color="#ffff00" bgColor="#010101" fontSize="10">
    <text lang="eng">
      <p bold="true">&lt;= leading fnct is '1'</p>
    </text>
  </region>
  <region id="12" left="216" top="1449" width="130" height="28" color="#000000" bgColor="#ffffff" fontSize="10">
    <math>
      <input>
        <e type="operand">n</e>
        <e type="operand">x</e>
        <e type="operand">β</e>
        <e type="function" args="2">p</e>
        <e type="function" preserve="true" args="1">rows</e>
        <e type="operator" args="2">≡</e>
      </input>
    </math>
  </region>
  <region id="13" left="216" top="1476" width="43" height="24" color="#000000" bgColor="#ffffff" fontSize="10">
    <math>
      <input>
        <e type="operand">n</e>
        <e type="operand">3</e>
        <e type="operator" args="2">:</e>
      </input>
    </math>
  </region>
  <region id="14" left="36" top="1521" width="192" height="66" color="#000000" bgColor="#e1ffff" fontSize="10">
    <math fractionType="none" decimalPlaces="6">
      <description active="true" position="Top" lang="eng">
        <p>2. Collect  the model function</p>
      </description>
      <input>
        <e type="operand">x</e>
        <e type="operand">β</e>
        <e type="function" args="2">f</e>
        <e type="operand">x</e>
        <e type="operand">β</e>
        <e type="function" args="2">p</e>
        <e type="operand">n</e>
        <e type="function" preserve="true" args="2">el</e>
        <e type="operand">n</e>
        <e type="operand">1</e>
        <e type="operand">3</e>
        <e type="function" preserve="true" args="4">sum</e>
        <e type="operand">1</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="3">line</e>
        <e type="operator" args="2">:</e>
      </input>
    </math>
  </region>
  <region id="15" left="279" top="1548" width="253" height="33" color="#000000" bgColor="#ffffff" fontSize="10">
    <math>
      <description active="true" position="Top" lang="eng">
        <p>view the expanded linfit model</p>
      </description>
      <input>
        <e type="operand">x</e>
        <e type="operand">β</e>
        <e type="function" args="2">f</e>
      </input>
      <result action="symbolic">
        <e type="operand">β</e>
        <e type="operand">2</e>
        <e type="function" preserve="true" args="2">el</e>
        <e type="operand">x</e>
        <e type="function" preserve="true" args="1">ln</e>
        <e type="operator" args="2">*</e>
        <e type="operand">β</e>
        <e type="operand">3</e>
        <e type="function" preserve="true" args="2">el</e>
        <e type="operand">x</e>
        <e type="function" preserve="true" args="1">sqrt</e>
        <e type="operator" args="2">*</e>
        <e type="operator" args="2">+</e>
        <e type="operand">β</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="2">el</e>
        <e type="operator" args="2">+</e>
      </result>
    </math>
  </region>
  <region id="16" left="36" top="1629" width="418" height="31" color="#000000" bgColor="#80ffff" fontSize="14">
    <text lang="eng">
      <p bold="true">Part 2: Execute the "linfit" solver</p>
    </text>
  </region>
  <region id="17" left="72" top="1674" width="97" height="36" color="#000000" bgColor="#ffffff" fontSize="10">
    <math>
      <description active="true" position="Right" lang="eng">
        <p>initialise the solver</p>
      </description>
      <input>
        <e type="operand">β</e>
        <e type="operand">1</e>
        <e type="operand">1</e>
        <e type="operand">1</e>
        <e type="operand">1</e>
        <e type="operand">3</e>
        <e type="function" preserve="true" args="5">mat</e>
        <e type="function" preserve="true" args="1">transpose</e>
        <e type="operator" args="2">:</e>
      </input>
    </math>
  </region>
  <region id="18" left="72" top="1728" width="179" height="26" color="#000000" bgColor="#ebd7ff" fontSize="10">
    <math decimalPlaces="6">
      <description active="true" position="Right" lang="eng">
        <p>&lt;= assign/create the vector gradient</p>
      </description>
      <input>
        <e type="operand">x</e>
        <e type="operand">β</e>
        <e type="function" args="2">φ</e>
        <e type="operand">x</e>
        <e type="operand">β</e>
        <e type="operand">f</e>
        <e type="operand">n</e>
        <e type="function" args="4">φ</e>
        <e type="operator" args="2">:</e>
      </input>
    </math>
  </region>
  <region id="19" left="72" top="1764" width="193" height="26" color="#000000" bgColor="#e1ffff" fontSize="10">
    <math fractionType="fraction" decimalPlaces="4">
      <description active="true" position="Right" lang="eng">
        <p>assign the solver</p>
      </description>
      <input>
        <e type="operand">β</e>
        <e type="operand">XY</e>
        <e type="operand">f</e>
        <e type="operand">φ</e>
        <e type="operand">β</e>
        <e type="function" args="4">βlinfit</e>
        <e type="operator" args="2">:</e>
      </input>
    </math>
  </region>
  <region id="20" left="72" top="1800" width="105" height="63" color="#000000" bgColor="#ffffff" fontSize="10">
    <math>
      <description active="true" position="Right" lang="eng">
        <p>coefficients of the model vector fnct's</p>
      </description>
      <input>
        <e type="operand">β</e>
      </input>
      <result action="numeric">
        <e type="operand">4.183</e>
        <e type="operand">0.996</e>
        <e type="operand">0.079</e>
        <e type="operator" args="1">-</e>
        <e type="operand">3</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="5">mat</e>
      </result>
    </math>
  </region>
  <region id="21" left="72" top="1872" width="157" height="67" color="#ffff00" bgColor="#010101" fontSize="10">
    <math>
      <description active="true" position="Right" lang="eng">
        <p>recast the solved "linfit model"</p>
      </description>
      <input>
        <e type="operand">x</e>
        <e type="operand">β</e>
        <e type="function" args="2">fit</e>
        <e type="operand">1</e>
        <e type="operand">x</e>
        <e type="function" preserve="true" args="1">ln</e>
        <e type="operand">x</e>
        <e type="function" preserve="true" args="1">sqrt</e>
        <e type="operand">3</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="5">mat</e>
        <e type="operand">β</e>
        <e type="operator" args="2">*</e>
        <e type="operator" args="2">:</e>
      </input>
    </math>
  </region>
  <region id="22" left="72" top="1944" width="240" height="158" color="#000000" bgColor="#ffffff" fontSize="10">
    <plot type="2d" render="lines" scale_x="1.4641" scale_y="0.59049" scale_z="0.864536409000001" rotate_x="0" rotate_y="0" rotate_z="0" transpose_x="-83" transpose_y="-47" transpose_z="0">
      <input>
        <e type="operand">data</e>
        <e type="operand">x</e>
        <e type="operand">β</e>
        <e type="function" args="2">fit</e>
        <e type="operand">2</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="4">sys</e>
      </input>
    </plot>
  </region>
</regions>